System and method for registration of cubic fisheye hemispherical images

ABSTRACT

Systems and methods for registering hemispheric images obtained using a fisheye lens for panoramic viewing, relating to spatial alignment and color balancing of complement half-side images derived from the hemispheric images. The systems and methods for spatial alignment determine the displacement of the digitized recorded images using a single translation and rotation model of distortion, which further evaluates the center of projection, and the extent of translation and rotation. The system and method for color balancing iteratively increase (or decrease) the values of pixels near the edge (discontinuity) between the two half-sides of an image, each time taking into consideration the average difference of pixel color across a smaller strip of the edge near the pixel. This invention removes the color differential but does not remove any detail underlying the features of the image.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation of U.S. application Ser. No. 10/433,906, filed Oct. 15, 2003 (now U.S. Pat. No. 7,412,091, issued Aug. 12, 2008), which is a national stage application of International Application No. PCT/CA01/01755, filed Dec. 7, 2001, which claims the benefit of U.S. Provisional Application No. 60/279,662, filed Mar. 30, 2001.

FIELD OF THE INVENTION

This invention relates to systems and methods for registering hemispheric images for panoramic viewing. In particular it relates to spatial alignment and colour balancing of images.

BACKGROUND AND GENERAL DESCRIPTION

Fisheye lens has experienced success in a number of applications involving panoramic or wide field-of-view applications. Such would include cinematography (U.S. Pat. No. 4,070,098), motionless surveillance (U.S. Pat. No. RE036,207), and image-based virtual reality (U.S. Pat. No. 5,960,108).

The advantage of fisheye projection is its large field of view compared to conventional rectilinear film. Images of field of view up to 220 degrees and beyond may be obtained with a fisheye lens. It has been speculated that a field of view infinitesimally less than 360 degrees is also obtainable, although the precise practical application of a lens of this type may be limited. In contrast, a conventional camera would require a rectilinear image recording surface of infinite dimensions for even 180 degrees of field of view.

As long ago as 1964, formal studies have been made of the optical characteristics of the fisheye lens (Kenro Miyamoto, “Fish eye lens”, Journal of Optical Society of America, 54:1060-1061, 1964). In 1983, Ned Greene suggested the use of fisheye images to generate environmental maps (Ned Greene, “A Method for Modeling Sky for Computer Animation”, Proc. First Int'l Conf. Engineering and Computer Graphics, pp.297-300, 1984).

In 1986, the use of perspective mapping for a fisheye image, projecting the latter into the sides of a rectangular box as an environmental map was introduced by Greene (Ned Greene, “Environmental Mapping and Other Applications of World Projections”, IEEE Computer Graphics and Applications, November 1986, vol. 6, no. 11, pp. 21-29). Greene took a 180 degree fisheye image (a fisheye environmental map) and projected it onto the six sides of a cube for perspective viewing.

Producing high-quality panoramic imaging using Greene's approach poses a number of difficulties. Each hemispheric image produces four half-sides of a cube, in addition to a full side, which require registration with its complement from the other hemispheric image. Registration of the half-images has two associated problems: spatial alignment and colour balancing.

Where the lens has a field of view greater than 180 degrees, the corresponding half-sides (as de-warped from the raw source image) require spatial alignment due to possible rotational and translational distortion. Furthermore, the image recording device would not necessarily capture the images precisely in the same area on the recording surface. Conventional methods of aligning the sides of such images are essentially manual in nature; even if assisted by graphics software, the process requires human intervention to line up the edges of half-sides based on the features in the images. This is often imprecise and difficult due to the multiple sources of distortion. There is a need for a more automatic process which reduces human intervention.

Chromatically, each half-side must be aligned relative to its complement. For example, the recording device may have the same exposure times and aperture openings for each of the two images despite the fact that the amount of light recorded for each image differs, as where the lighting changed between the two capture events, or alternatively, where the exposure times are different for each image despite equivalent lighting conditions. If aperture size is controlled automatically, further mismatches may result. As a result the complementary edges of half-images have different colour intensity but generally the same detail and level of noise.

Existing methodology for colour balancing tend to average the colour values of pixels in the relevant neighbourhood of the transition. This has the desired effect of bringing the difference across the transition into line. However, the disadvantage concerns the concomitant loss of detail. As a result there is a perceptible blurring across the region of the transition. When the width of the overlapped region is narrow as compared to lighting imbalance, the transition may appear too abrupt.

SUMMARY OF THE INVENTION

This invention provides systems and methods for registering hemispheric images obtained using a fisheye lens for panoramic viewing. In particular it relates to spatial alignment and colour balancing of complement half-side images derived from the hemispheric images.

The system and method for spatial alignment determines the displacement of the digitized recorded images using a single translation and rotation model of distortion, which further evaluates the centre of projection, the distance in pixel for 180 degree in the recorded image, the extent of translation and rotation.

The system and method for colour balancing iteratively increase (or decrease) the values of pixels near the edge (discontinuity) between the two half-sides of an image, each time taking into consideration the average difference of pixel colour across the edge of a strip of decreasing width. This invention removes the colour differential but does not remove any detail underlying the features of the image.

BRIEF DESCRIPTION OF THE DRAWINGS

Embodiments of the invention will be described by way of example and with reference to the drawings in which:

FIG. 1: A diagram showing the relationship between the angle subtended by a point and the radial position on the recorded source image.

FIG. 2: A diagram showing the mapping from the hemispheric image to 5 sides of the cube.

FIG. 3: A diagram indicating how the Cartesian coordinates of a point correspond to its spherical coordinates.

FIG. 4: A diagram of the two sides of an edge (discontinuity in colour regimen) as represented by two 2-dimensional arrays.

DETAILED DESCRIPTION OF THE INVENTION

The description given below assumes the use of an ideal fisheye lens, meaning that the lens possesses negligible if no radial distortion: the radial position recorded in the source image for a point imaged is proportional to the zenith angle the point subtends with respect to the axis of projection (see FIG. 1). However, a person knowledgeable in the art would be able to extend the invention disclosed below to the case where there is radial distortion. Adjustment to take into consideration lens characteristics are clear once the characteristics are determined.

Furthermore, the fisheye lens has a field of view greater than 180 degrees.

Overview

In the present invention, an environmental map of cube consisting of 6 canonical views is reconstructed from two hemispherical fisheye images, as Greene had suggested, each with field of view equal or greater than 180 degrees. Each fisheye image is manipulated to generate 1 complete side and 4 half-sides of a cube; and the half sides are then integrated to form 4 complete sides; all six sides as the environmental map are then placed together for viewing as a complete cube texture map using computer graphic techniques. See FIG. 2 for an overall depiction of how a hemispherical image may be decomposed into the various sides of a cube.

The spherical coordinates, ρ, θ, φ, expressed in Cartesian coordinates, x, y, z, of a point in space, are as follows (see FIG. 3):

$\rho = \sqrt{x^{2} + y^{2}}$ $\theta = {\tan^{- 1}\left( \frac{y}{x} \right)}$ $\phi = {\tan^{- 1}\left( \frac{\rho}{z} \right)}$ θ is the azimuthal angle, and φ, the zenith angle. The effect of the fisheye lens is such that each point is projected onto a hemisphere of diameter D, centred at the origin, with the projective axis being the z-axis. The point on the hemisphere when subsequently projected onto the x-y plane, corresponds to the source image coordinate (u, v) as captured by the image recording device. As a result, the value of the pixel at location (u, v) of the source image is assigned as the value of the pixel on the environmental map that corresponds to the (x, y, z) coordinate on the cube.

Given the above, for any point on the hemi-cube, the coordinate (u, v) in the source image (the environmental map) can be evaluated from θ and φ, corresponding to the planar coordinates of the point of intersection of a ray extending from the origin to the hemi-cube with the hemi-sphere, as follows.

$u = {\frac{\phi\; D}{\pi}{\cos(\theta)}}$ $v = {\frac{\phi\; D}{\pi}{\sin(\theta)}}$

Since the lens is ideal, the radial distance from the origin of any point (u, v) on the x-y plane should be linear in the zenith angle φ. (Take the square root of the squares of u and v, which results in a linear function of φ.) Given a digitized hemispheric fisheye source image formed as described above, the starting point is to determine the centre of projection (u_(c), v_(c)) on the source image array, and the length in image pixel, D, equal to a view of 180 degrees in zenith angle. In practice, the centre of projection, (i.e. the point where φ=0°) is not always the centre of the source image array. One would then determine the image array corresponding to each of the sides by determining for each point in the corresponding image array the (x, y, z) coordinate, compute the corresponding (u, v) coordinate on the source image, and estimate the value of the pixel at (u, v) using some sort of interpolation scheme based on the neighboring pixel values. In the expression for (u, v), the coordinate must be translated by the centre of projection (u_(c), v_(c)) in the source image.

At this point the images require colour balancing in order to remove visible discrepancies in colour between any two complementary half-images along the separating edge.

With the geometry aligned and colour balanced, the visible interface between the hemi-cubes disappears, resulting in a complete cube texture map comprising the six full-sides. The stacked face arrangement and ordering is preferably chosen to comply with a VRML cube primitive, for possible display on the world wide web.

Standard computer graphic techniques may be used to render the six-sided cube centred at the origin. The viewpoint is fixed at the origin, and the view direction along the positive z axis. Using a pointing device such as a mouse, a viewer can vary the yaw and pitch of the view by rotating the cube around the Y-followed by a rotation about the X-axis. The faces are mapped to the cube using perspective texture mapping, and as the user switches scenes, the texture map changes.

Preferred embodiments of this invention use one or more electronic processors with graphics display devices and other human interface devices. The processors may be a native device or resident on an electronic network such as a local area network, an intranet, an extranet, or the Internet. The source images may be derived from a digital camera, such as those using CCD cells, or from digitized images of photographs taken with a conventional camera outfitted with a fisheye lens, amongst others.

Geometric Alignment

Spatial, or geometric, alignment of the two hemispheric images proceeds on the basis that there is an overlapping annular region between the two images. Given a lens with a half field of view of ν=105 degrees, the annulus of overlap would be 30 degrees (two times 15 degrees on each side) if the hemispheric images were taken exactly 180 degrees apart. In this case, each annular region beyond the ring denoting 180 degrees would be a mirror image of the annular region of the same dimension in the other hemispheric image (reflected across the 180° ring). Thus in FIG. 1 the region consisting of the two outermost concentric annuli in each of the mapping of the two fisheye lens is a mirror image of the another.

In practice, the camera may experience displacement from the ideal position for taking the reverse image after recording the initial one: rotation about each of the 3 Cartesian axes and translational movement in the x-y plane may occur. The model of distortion assumed by this invention however postulates 2 types of distortion: a translation in the x-y plane and a rotation about the projection z-axis. The image recorders are assumed to be directed precisely 180 degrees apart with coincidental projection axes, the z-axis. The source image planes are assumed to be coplanar, but the image as recorded experiences a possible translational and rotational displacement.

The approach disclosed by this invention is to align 3 points in the overlapping region of one hemispherical image with the same 3 points in the other image using the above mentioned model of distortion. This method aligns the three pairs of points such that the zenith angle φ from the centre of projection is the same for any two pairs of points. Once aligned, the centre of projection and the distance in pixel for 180 degrees field of view may be computed. Knowledge of these two allows the cubic faces of the environmental map to be generated.

The following defines the notation adopted:

-   -   3 non-collinear points p₁(x₁, y₁), p₃(x₃, y₃), and p₅(x₅, y₅) in         the overlapping annular region of a hemispheric source image T₁         with the centre of projection at q₁ and the complement set of         points p₂(x₂, Y₂), P₄(₄, y₄), and p₆(x₆, y₆) with its centre of         projection at q₂ for the reverse image T₂ (Since T2 is the         reverse image, it becomes necessary to reflect T₂ about the         vertical axis defined by the azimuthal angle φ being 90 or 180         degrees, to yield T₂′. For the purpose of simplicity though,         refer henceforth to T₂′ as T₂;     -   the angle that each source image is rotated about its axis of         projection relative to the earth's horizon: α₁ and α₂, for T₁         and T₂ respectively. For the purpose of simplicity, assume that         α₂ is zero, that is to say, the second image is not rotated         relative to the horizon;     -   the true center of projection: Q;     -   the diameter in pixels equivalent to 180 degrees field of view:         D.

Given that the images T₁ and T₂ are precisely 180 degrees apart, any point p₁ will form the same azimuthal angle θ relative to q₁, the centre of projection in its own image, as its complement P₂ would with q₂. Therefore, if the two centres of projection were the same point Q, then the z component of the vector cross product of (p₁−Q) and (P₂−Q) would equal zero since these two vectors are parallel.

The centre Q may be found assuming that displacement consists solely of a translation Δx and Δy, and a rotation of angle α. One embodiment of this approach assumes that q₂ is fixed as the centre Q, and determines the displacement (Δx, Δy, α) of q₁, resulting in q₁′, such that the z components of the 3 vector cross products of (q₁′−p₁) and (q₂−p₂), (q₁′−p₃) and (q₂−p₄), and (q₁′−p₅) and (q₂−p₆) are all zero. (The cross products are all zero vectors since the initial vectors have zero z component resulting in no planar components to the resultant cross products.)

This essentially sets up a system of 3 equations (the z components of the cross products) in 3 unknowns (Δx, Δy, α), which a person skilled in the art is able to solve. For example, the Newton-Raphson iterative method for finding roots of systems may be used. An initial approximation to the solution for (Δx, Δy, α) may be the translational difference between p₁ and P₂, and the angle zero.

The result of this computation is the 3-tuple, (Δx, Δy, α). If p_(i)′ is the result of translating and rotating p_(i), each complementary pair of vectors (q₂−p_(i)′) and (q₂−p_(j)), where i=1, 3, and 5, and j=2, 4, and 6 respectively, are parallel, thus the lines joining each complement pair of p_(i)′ and p_(j) must intersect at the same point q₂ (same position as Q and q₁′). As a result, the centres of projection in each hemispheric image (referred to as (u_(c), v_(c)) for each hemispheric image earlier), q₁ and q₂, may be established.

Furthermore, it is now possible to determine D, the number of pixels equal to 180 degrees field of view, which is used for evaluating the value of coordinates (u, v) in the environmental map. Since the sum of the zenith angle φ of a point in one image added to the zenith angle of the same point in the other hemispheric image is 180 degrees, the zenith angle 180 corresponds to the average of the radial distance for any two complement points. Therefore this allows D to be computed.

Given the centre of projection, (u_(c), v_(c)), for each hemispheric image, and the distance in pixel D for 180 degrees, the 6 perspective sides of the cubic environmental map may now be generated.

In one preferred embodiment of this invention, the user selects for display annular regions from both hemispheric images. Each annulus covers typically the overlapping region between the two source images, subtended by the zenith angles between 270−ν and 90+ν, where ν is the half-field of view. For example, if ν is 105 degrees (i.e. the full field of view is 210 degrees), the annulus region corresponds to the field of view between 165 and 195 degrees. Once the two annuli of interest are displayed, typically at the same time, the user may then select a common point in each annulus. This invention also includes the variation where the user merely has to denote the approximate area where the second of the pair of points is located in the second image. The electronic processor would utilize imaging processing techniques such as a convolution mask to locate the precise point in the second image corresponding to the first point in the first image.

This invention further includes variations where automatic determination of the 3 common points are performed without human intervention.

Although it is theoretically possible to select all three points in the same nearby neighbourhood, the preferred approach is to select points which are widely spread apart, typically in different thirds of the annular region. Therefore, display of only a common sector of the annuli typically occurs for the point selection process. Displaying smaller areas permits greater spatial resolution and more precise location within the areas.

Once the common points are identified, the electronic processor carries out computation of the centre of projection (u_(c), v_(c)) for each hemispheric image and the distance D in pixel for 180 degrees field of view. The environmental map of the cubic sides are then ready to be fully generated given these parameters and the rotation α.

Colour Balancing

In addition to removing spatial distortions as mentioned above, two complementary half-images as combined typically require colour balancing due to different prevailing lighting conditions or imaging recorder characteristics. This invention further discloses methods and apparatus for colour balancing across a discontinuity (or edge) in an image that does not remove detail from the image as a whole.

To illustrate the method by example, let one side of the edge be represented by a 2 dimensional array A_(j,i), and the other side B_(j,i), both of dimensions m+1 rows (indexed 0 to m) and n+1 columns (indexed 0 to n). The 0-th row of both arrays lie adjacent to the edge as shown in FIG. 4, with rows further away from the edge in increasing index of row. Although only one value for each pixel is indicated here, this method will work for any color model including RGB, CMY, YIQ, HSV, and HLS. These colour models are well known in the field of computer graphics. For example, in the case of an RGB colour model, there would be 3 arrays, ₁A_(j,i), ₂A_(j,i), and ₃A_(j,i), corresponding to each of the red, green, and blue colour values; the same would apply for the B arrays.

This invention makes use of the concept of “imbalance waves” of varying frequency which runs additively parallel and adjacent to the edge for a certain depth normal to the edge. An iterative approach is adopted whereby a different filter is applied at each step to the edge region. The analogy is that each filter is a low pass filter, tuned to a distinct “frequency” level, which eliminates imbalance waves of frequency at that level and greater. During each iteration, a filter of a lesser frequency is used until a frequency of 1 is reached.

Generally, at each iteration, a filter of level w=2^(h) is applied, where h decreases from m to 0. The value w is a measure of the width of a strip on the edge, the average difference, in colour value across the edge on the strip, is permitted to influence colour values during the steps constituting the particular iteration. The value m is chosen such that 2^(m+1) is less than the width of a full side of the cube, which will become clear later in this discussion. Instead of powers of two, which has the advantage of being intuitive and can enable m to decrease to the value 0 quickly as the logarithm of w, other schemes for decreasing w will also suffice.

During an iteration, the difference at column i across the edge between any two adjacent A_(0,i) and B_(0,i) is evaluated as Δ_(i.) Δ_(i) =A _(0,i) −B _(0,i) i=0, 1, . . . , n

In cases where there is high contrast in difference across an edge, a slightly different method for computing Δ_(i) is helpful to avoid “bleeding” of the high intensity value of the one side into the lower intensity opposite side in the balanced image. This approach assumes the images as slightly overlapping in the rows and slightly offset in columns (which is typically the case in reality), and uses both images. If it is assumed that there are x overlapping rows, one way of computing the value of Δ_(l) is to use a small region across the boundary in one image and centred on A_(0,i) (or B_(0,i), no difference whichever is used) as a square mask (typically 3 elements by 3 elements). The region is overlaid on a series of corresponding areas in the other image, the other areas within a typically square domain centred on the ′A_(0,i) in the other image, where ′A_(0,i) is the corresponding pixel for A_(0,i) in the original image. Each side of the domain has size typically twice the overlapping width (plus one). The sums of the squares of the differences between the values of the mask and the underlying pixels are calculated as the mask is translated across the surface of the domain. One embodiment calculates Δ_(l) as the average of the arithmetic difference of the mask the underlying pixels of the image centred at the place of the minimal sum of squares.

If the overlap is greater than two, then the mask may be chosen be larger than 3×3, up to the point where it may be limited in translational freedom of only one pixel in any of the 8 directions.

An average of the differences, Γ_(i), for any column i on the edge is taken for 2w+1 differences centred on i, with wraparound values for any difference with index exceeding the bounds of the difference array.

$\Gamma_{i} = \frac{\sum\limits_{k = {i - w}}^{i + w}\Delta_{k}}{{2\; w} + 1}$ where k=k+(n−1), if k<0, and k=k−(n−1), if k>n, for the purpose of the summation. Clearly, this is just one way of calculating the average of the difference, and the invention is not restricted to this one way, although this results in computational efficiency. Other variations could weight the differences by a window appropriate to the type of image, e.g. a Gaussian window centred on i. Alternatively, there may be a threshold for discounting those differences which exceed the threshold. For example, the presence of a colour burst may be one instance where a precipitous drop in value across an edge should be discounted. Visually, such a measure would have little impact on the resultant reconstructed image.

At this point, the value of the side with the lower value is adjusted upwards by a variation term which is a function of the average difference across the edge. The adjustment used decreases with distance from the edge. A preferred embodiment then redefines all the elements of the arrays, A_(j,i) and B_(j,i), as follows:

$A_{j,i} = {A_{j,i} - {\frac{1}{2}\Gamma_{i}\frac{w - j}{w}}}$ $B_{j,i} = {B_{j,i} + {\frac{1}{2}\Gamma_{i}\frac{w - j}{w}}}$

Although the above weights the average difference decreasingly with distance on a linear basis from the edge (w−j), this invention is not limited to this weighting scheme. One variation is where a faster drop-off, e.g. an inverse exponential, is required by the character of the image.

The way that Δ_(i) is calculated influences the effect of a high contrast edge. Without seeking out the minimum sum of square of difference, a high contrast edge will change the value of both sides of the edge more significantly than if the difference where the minimum sum of square occurred was used.

Typically, the iteration ends when w reaches the value of one (filter of frequency 1). However, it is possible to impose a threshold beyond which no further iterations are carried out. One approach is based on the average difference across the edge. Once the average difference drops below a threshold, the iterative steps end. In this case the threshold may be taken relative to a measure of the noise in the neighbourhood of the edge.

Preferred embodiments of this invention also involve computer program products which are recorded on a machine-readable media. Such computer program products when loaded into an apparatus that can read and execute the instructions borne by the program products in order to carry out the functions of geometric alignment or colour balancing denoted above.

The system, methods, and products above may operate modularly or in serial. For example, the same system may perform both alignment and colour balancing, or two separate systems operate serially on the digitized images. These may be used in standalone configurations or in a networked environment cooperating with many other applications, systems, databases etc. for capturing and processing. Environmental maps may be generated from colour balanced images.

It will be appreciated that the description above relates to the preferred embodiments by way of example only. Many variations on the device, method, and computer program product for delivering the invention will be understood to those knowledgeable in the field, and such variations are within the scope of the invention as described, whether or not expressly described. 

The invention claimed is:
 1. A non-transitory computer readable medium having processor-executable instructions configured to cause one or more processors to perform the steps of: (A) receiving a plurality of images in at least one input file, an edge being defined where a first given row of a first of said images meets a second given row of a second of said images and said edge is normal to a plurality of columns of said images; (B) calculating a width of a strip based on an iteration value, said strip being oriented parallel to said edge; (C) calculating a plurality of difference values across said edge, one of said difference values for each respective one of said columns, each respective one of said difference values being based on a plurality of first color values of said first image inside said strip and a plurality of second color values of said second image inside said strip with said strip centered about said respective column; (D) adjusting said first color values and said second color values in each of said columns by a weighted scheme based on both (i) said respective difference value and (ii) a respective distance normal from said edge; (E) repeating the steps (B) through (D) for each of a plurality of said iteration values, each successive one of said iteration values decreasing said width of said strip; and (F) storing said images as adjusted in an output file.
 2. A non-transitory computer readable medium according to claim 1, wherein step (B) comprises the sub-step of: calculating said width as w=2^(h) where h decreases from m to 0, and where m is chosen such that 2^(m+1) is less than a width of a full side of a viewing cube.
 3. A non-transitory computer readable medium according to claim 2, wherein step (C) comprises the sub-step of: calculating said difference values as Δ_(i)=A(0,i)−B(0,i) over i=0, . . . , n where Δ_(i) are the difference values at column i across the edge, where one side of the edge is represented by two dimensional array A(j,i), and wherein the other side of the edge is represented by two dimensional array B(j,i).
 4. A non-transitory computer readable medium according to claim 3, wherein step (C) comprises the sub-step of: calculating said respective difference value as a sum of a plurality of squares of a plurality of differences between (i) a plurality of mask values of a mask and (ii) said first color values and said second color values aligned with said mask.
 5. A non-transitory computer readable medium according to claim 4, wherein step (D) comprises the sub-step of: calculating an average of said difference values in said strip as: $\Gamma_{i} = \frac{\sum\limits_{k = {i - w}}^{i + w}\Delta_{k}}{{2w} + 1}$ where k=k+(n−1), if k<0, and k=k−(n−1), if k>n, for the purpose of the summation.
 6. A non-transitory computer readable medium according to claim 5, wherein step (D) further comprises the sub-step of: adjusting said first color values by: ${A_{j,i} = {A_{j,i} - {\frac{1}{2}\Gamma_{i}\frac{w - j}{w}}}};$ and adjusting said second color values by: $B_{j,i} = {B_{j,i} + {\frac{1}{2}\;\Gamma_{i}\frac{w - j}{w}}}$ where A_(j,i) has a higher value than B_(j,i).
 7. A non-transitory computer readable medium according to claim 5, further comprising the step of: ending said repeating in response to said average of said difference value becoming less than a threshold.
 8. A non-transitory computer readable medium according to claim 7, wherein said threshold comprises a measure of noise in a neighborhood of said edge.
 9. A non-transitory computer readable medium according to claim 1, wherein said processor executable instructions comprise a stand-alone program.
 10. A non-transitory computer readable medium according to claim 1, wherein said processor executable instructions comprise comprises an executable module within a program.
 11. A non-transitory computer readable medium according to claim 1, wherein said edge is formed along a straight line.
 12. One or more processors including structure configured to: (A) receive a plurality of images in at least one input file, an edge being defined where a first given row of a first of said images meets a second given row of a second of said images and said edge is normal to a plurality of columns of said images; (B) calculate a width of a strip based on an iteration value, said strip being oriented parallel to said edge; (C) calculate a plurality of difference values across said edge, one of said difference values for each respective one of said columns, each respective one of said difference values being based on a plurality of first color values of said first image inside said strip and a plurality of second color values of said second image inside said strip with said strip centered about said respective column; (D) adjust said first color values and said second color values in each of said columns by a weighted scheme based on both (i) said respective difference value and (ii) a respective distance normal from said edge; (E) repeat the steps (B) through (D) for each of a plurality of said iteration values, each successive one of said iteration values decreasing said width of said strip; and (F) store said images as adjusted in an output file.
 13. One or more processors according to claim 12, wherein the one or more processors includes structure configured to, in calculating the width of the strip: calculate said width as w=2^(h) where h decreases from m to 0, and where m is chosen such that 2^(m+1) is less than a width of a full side of a viewing cube.
 14. One or more processors according to claim 13, wherein the one or more processors includes structure configured to, in calculating the plurality of difference values: calculate said difference values as Δ_(i)=A(0,i)−B(0,i) over i=0, . . . , n where Δ_(i) are the difference values at column i across the edge, where one side of the edge is represented by two dimensional array A(j,i), and wherein the other side of the edge is represented by two dimensional array B(j,i).
 15. One or more processors according to claim 14, wherein the one or more processors includes structure configured to, in calculating the plurality of difference values: calculate said respective difference value as a sum of a plurality of squares of a plurality of differences between (i) a plurality of mask values of a mask and (ii) said first color values and said second color values aligned with said mask.
 16. One or more processors according to claim 15, wherein the one or more processors includes structure configured to, in adjusting said first color values and said second color values: calculate an average of said difference values in said strip as: $\Gamma_{i} = \frac{\sum\limits_{k = {i - w}}^{i + w}\Delta_{k}}{{2w} + 1}$ where k=k+(n−1), if k<0, and k=k−(n−1), if k>n, for the purpose of the summation.
 17. One or more processors according to claim 16, wherein the one or more processors includes structure configured to, in adjusting said first color values and said second color values: adjust said first color values by: ${A_{j,i} = {A_{j,i} - {\frac{1}{2}\Gamma_{i}\frac{w - j}{w}}}};$ and adjusting said second color values by: $B_{j,i} = {B_{j,i} + {\frac{1}{2}\;\Gamma_{i}\frac{w - j}{w}}}$ where A_(j,i) has a higher value than B_(j,i).
 18. One or more processors according to claim 16, wherein the one or more processors includes structure configured to: end said repeating in response to said average of said difference value becoming less than a threshold.
 19. One or more processors according to claim 18, wherein said threshold comprises a measure of noise in a neighborhood of said edge.
 20. One or more processors according to claim 12,wherein said edge is formed along a straight line. 